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Abstract An approximate solution of the position-dependent mass Dirac equation with the 
Hulthen potential is obtained in D-dimensions within frame work of an exponential approx- 
imation of the centrifugal term. The relativistic energy spectrum is worked out using direct 
transformation method; the two-component spinor wavefunctions are obtained in terms of the 
Jacobi polynomials. Dependence of the energy levels on some parameters is discussed. The 
results obtained are in good agreement with previous works. 
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1 Introduction 

The relativistic Dirac equation plays an important role in the understanding of many quantum 
mechanical problems. The equation which describe the motion of a spin-^ particle has been 
solved with various physical potentials. For instance, the solutions of the Dirac equation 
with the Coulomb potential have been obtained in ID [1, 2], 2D [3] and 3-D [4-6]. Also, the 
qualitative properties of the Dirac particle in a central potential and the discrete eigenvalues of 
the radial Dirac operator have been discussed [7, 8]. However, with the intreest in the higher 
dimensional field theory, the multidimensional relativistic and non-relaivistic equations have 
been studied by many authors. To mention a few, the D-dimensional Schrodinger has been 
studied with the Coulomb potential [9], pseudoharmonic potential [10], Hulthen potential [11] 
and Poschl- Teller potential [12]. Also, various potentials have been studied with both the 
D-dimensional klein-Gordon [13-16] and Dirac [17-20] equations. 

Moreover, the Hulthen potential is one of the important short-range potentials which 
behaves like a Coulomb potential for small values of r and decreases exponentially for large 
values of r. The Hulthen potential has received extensive study in both relativistic and 
non-relativistic quantum mechanics [11, 13-16, 21, 22]. Unfortunately, quantum mechanical 
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equations with the Hulthen potential can be solved analytically only for the s-states [23-25]. 
However, recent study [26] has used an exponential approximation for the centrifugal term in 
order to obtain the solution of the I 7^ states of the Hulthen potential. This approximation 
has been empolyed to obtain the solutions of both the relativistic [14] and non-relativistic [11] 
Hulthen potential in D-dimensions. On the other hand, the used of position-dependent mass 
has produced interesting results in many bound-state problems. This approach has been used 
by many authors which include [27-30]. The present study is of two- fold. The first purpose 
is to study the bound state of the Dirac-Hulthen problem in D-dimensions. Secondly, we to 
investigate if the dimensional degeneracy observed in the non-relativistic Hulthenspectrum 
[11] also occur in the case of the relativistic equations. 

The paper is organized as follows. In section 2, we present the Dirac-Hulthen problem in 
D-dimensions. The solutions to the two-component wavefunctions and the relativistic energy 
spectrum are obtain in section 3. In section 4, the dimensional degeneracy for both Klein- 
Gordon and the Dirac equations are discussed and some concluding remarks are given in 
section 5. 



2 The Dirac equation in D-dimensions 

The D-dimensional Dirac equation with a central potential V(r) and position-dependent mass 
fi(r) can be written in natural units h = c = 1 as [18, 31] 

D 

HV(r) = E nrK y(r) where H = ^ UjPj + (3fi(r) + V(r) (1) 

i=i 

where E nK is the relativistic energy, {&j} and /3 are Dirac matrices, which satisfy anti- 
commutation relations 

djOfc + QfcQj = 25jkl 

ajf3 + /3aj = (2) 
«] = /3 2 = 1 

and 

Pi = -idj = -i^- l^j^D. (3) 

The orbital angular momentum operators Ljk, the spinor opertaors Sjj~ and the total angular 
momentum operators Jjk can be defined as follows: 

d . d 

Ljk = — Ljk = i x j~Q^ ^ Xfc 0x~' = ~^ k i = Jjk = Ljk + Sjk- 

D D D 

L 2 = Y,L%, S 2 = Y,S* k , J 2 = Y. J h> l^j<k^D. (4) 

j<k j<k j<k 

For a spherically symmetric potential, total angular momentum operator Jjk and the spin- 
orbit operator K = —/3(J 2 — 1? — S 2 + (D — l)/2) commutate with the Dirac Hamiltonian. 
For a given total angular momentum j, the eigenvalues of K are k = ±(j + (D — 2)/2); 
k = — (j + (D — 2)/2) for aligned spin j = I + \ and k = (j + (D — 2)/2) for unaligned spin 
j = I — \. Also, since V(r) is spherically symmetric, the symmetry group of the system is 
SO(D) group. 

Thus, we can introduce the hyperspherical coordinates [32-35] 
x\ = r cos 61 

x a = r sin^i . . . sin^Q,-! cos(/>, 2 ^ a ^ D — 1 (5) 
xd = r sin^i . . . sin6*D-2 sin0, 
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where the volume element of the configuration space is given as 

D D-l 

Y\ dxj = r D - l drdtt dn = Yl (smO j ) i - 1 dB j (6) 

3 = 1 J'=l 

with ^ r < oo, ^ 9k ^ 7r, k = 1, 2, . . . D — 2, 0^0^ 27r, such that the spinor 
wavefunctions can be classified according to the hyperradial quantum number n r and the 
spin-orbit quantum number k and can be written using the Pauli-Dirac representation 

( F nrK (r)Yf m (n D ) \ 
*n rK (r,n D ) = r"¥ (7) 
\ iG nrK (r)Yf m (Q D ) J 

where F nrK (r) and G 

n r K( r ) are the radial wave function of the upper- and the lower-spinor 
components respectively, Y-- m {Q,jj) and Yj m (Q£>) are the hyperspherical harmonic functions 
coupled with the total angular momentum j. The orbital and the pseudo-orbital angular 
momentum quantum numbers for spin symmetry I and and pseudospin symmetry t refer to 
the upper- and lower-component respectively. 

Substituting Eq. (7) into Eq. (1), and seperating the variables we obtain the following 
coupled radial Dirac equation for the spinor components: 



(j r + F nrK (r) = Kr) + E nrK - V(r)]G nrK (r) 
G nrK {r) = Mr) ~ E nrK + V(r)]F rirK (r) 



(8) 



(9) 



where k = ±(21 + D — l)/2. Further details of the derivation can be obtain from refs [36-38]. 
Using Eq. (8) as the upper component and substituting into Eq. (9), we obtain the follwoing 
second order differential equations 



d 2 k(k + 1) 
dr 2 r 2 

(f_ _ k(k-1) 
dr 2 r 2 



W) + E nrn ~ V(r)Mr) - E nrK + V{r) 



[H{r) + E nrK - V(r)Mr) - E nrK + V{r) 



dr 



dr 



W)+E nrK -V{r)\ 



Mr) , dV(r) \ ( 
dr ~r dr J \ 



A 

dr 



Hr)-E nrK + V(r)] 



F nrK (r) = 
(10) 

G nrK {r) = 
(11) 



We note that the energy eigenvalues in these equation depend on the angular momentum 
quantum number i and dimension D. However, to solve these equations, we shall use an 
approximation for the centrifugal barrier as discussed in the following section. 

3 Bound states of the Dirac-Hulthen problem in L>-dimension 

We start this section by defining the Hulthen potential as follows [11, 13, 14, 21, 22] 



V(r) 



-Za- 



1-e" 



(12) 



where a is the screening parameter and Z is a constant which is identified with the atomic 
number when the potential is used for atomic phenomenon. To solve Eq. (10), we first elimi- 



nate the last term by equating 



djMr) dV(r) 



dr 



-g^r- — 0, which gives the mass function 
Za 



V>(r) = no + 



1 - e- 



(13) 
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where /jLq is the intergation constant, and approximate the centrifugal term as follows [11, 26] 

1 



a 2 e~ ar 



r 2 (l-e- ar ) 2 ' 
Substituting Eqs. (12), (13) and (14) into Eq. (10) we have 



(14) 



d 2 K(n + l)a 2 e- ar 
dP 2 (1 - e~ ar ) 2 

Za 



Za 



1-e" 



E nrK + Za 



x A*o + 



E„ 



Za- 



■y _ e -ar '" rr ° ]_ _ e ~ar 

Taking the transformation s = e~ ar , Eq. (15) becomes 



1-e- 
F nrK {r) = 



d 2 Id k(k + 1) e 2 fh+fo ] p M _ n 

d S 2 + Sds S{l-S) 2 S 2 8(1-8)] ^ 



(15) 



(16) 



where 



^ 2 + ^a 2 -E 2 rK 2ZfiQ + Z 2 a 



a 



; Pi 



a 



and 02 



2ZE nrK + Z 2 a 



a 



If we seek the solution of the form 



then Eq. (16) becomes 



F nrK (s) = S £ (l " s) S U nrK ( S ), 

T + 2e 25 



(17) 



(18) 



(-(/3i + &) - (2e + 1)<5 - 5 2 + <5) (1 - s) + 5 2 - 5 - k(k + 1)' 



a(l -s) 2 



E^kOO = (19) 



One can choose S 2 — 5 — k(k + 1) = 0, i.e 8 = k + 1 (the positive root), such that Eq. (19) 
becomes 

s(l - s)U'* rK {s) + [(1 + 26) - (26 + 25 + 1»] E^s) 



^ + e - Ve 2 " (#i + ft)! [<f + e + V^ 2 1 (ft + fa)] U nrK (s) = 



(20) 



Eq. (20) is the well-known hypergeometric equation whose solution is given in form of the 
hypergeometric function [39] 



U nrK (s)= 2 F 1 5 + e- v / e 2_(/? 1+/ 3 2 ), s + e + y/e^-i^+lh); 1 + 2e; s 



(21) 



However, for large value of s, the solution in (20) diverges, thus preventing normalization. To 
avoid this, we set 



-n r = S + e - \A 2 -(/3i+/?2), n r = 0,1,2,.... 
from which, with the help of Eq. (17) we have the energy eigenvalues 

E nrK = -^L±M ± K+ f l + 1 V^o + a 2 A).-a 2 ( ?? + ^ 
2?7 277 



(22) 



(23) 
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where r\ = (n r + \n\ + l) 2 + Z 2 . To check the validity of the energy spectrum, we take the 
limit as a — >■ of Eq. (23) and this yields 



Ecol 



1 + 



(n r + \k\ + l) 2 



(24) 



which is the energy spectrum of the Coulomb-like potential in D-dimensions [17, 36]. Evi- 
dently, this follows from the fact that lim a ^o^( r ) = —y- 

Thus the unnormalized wavefunction for the upper-component can be written as 



F nrK (s) = C nrK s e (l - s) K+1 2 F 1 [-nr, k + e + 1 + v / e 2 -(/3i+/3 2 ); 1 + 2e; s 
Using the relation using the following definition of the Jacobi polynomial [39] 



T(n + a+ 1) 
n!r(l + a) 



F\ ( — n, a + b + n + 1; 1 + a; 



we arrive at 



F nrK {s) = C nrK s*(l - sT +l P^ 2K+1 )(1 - 2s) 



(25) 



(26) 



(27) 



where C HrK is the normalization constant. 

Moreover, substituting Eq. (27) into Eq. (8), one can easily obtain the lower-component as 

G nK (s) = A 1 (s)Pi 2 ;> 2K+1 )(1 - 2a) + A 2 (s)P^' 2k+2) (1 - 2s) (28) 

where 

B rir ^(l - s) K 



fi(r) + E nrK - V{r) 



(29) 



with B HrK a constant. 



4 Dependence of the relativistic Hulthen spectra on some 
parameters 

4.1 Case 1: Dependence of E n on a 

Recently, it was reported in [11] that the non-relativistic energy spectrum of the Hulthen 
potential intersect for some adjacent dimensions D(> 3), for some values of the screening 
parameter «j. Thus signifying a degeneracy (with respect to the dimension) in the energy 
spectrum of the non-relativistic Hulthen potential. In this section, we investigate if such 
degeneracy also occur in the relativistic energy spectums of the Hulthen potential. We start 
with the Klein-Gordon-Hulthen energy spectrum for equal vector and scalar potentials, which 
is given as [11, 13, 14] 

Ek„ = ± <!4p£vV-°<2^-7) (30) 

where n r is the radial quantum number, fj, is the mass and 

V! = (D + 2e-l)/2, 7 = 2Z/i-a(n r + ?;i) 2 and n = (n r + vxf + Z 2 . (31) 

Moreover, if we set Z = n = 1, and define a principal quantum number n = n r + i + 1, the 
energy eqution (30) takes the simple form (with the choice of a negative root) 

2ap x - 16 pi r— — 

^ = ^T7r-4T7r Vl6 " api - (32) 
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where p\ = (2n + D — 3) 2 . Similar to result of [11], the curve of Ek u against a (Figure 1) 
reveals the intersection of the energy levels for some adjacent dimensions. Thus indicating 
degeneracy (with respect to the dimensions) in the energy levels. 

We now consider the case of the Dirac-Hulthen spectrum. If we introduce a principal 
quantum number n = n r + |k| + 1, the energy equation (23) can be written (Z = p,® = 1) in 
the form 



**■ = - a( " 2te I T ^ - 2(^1) + ')(« + ')' - M» + 2) + 2]^ (33, 

where p 2 = (n + (D — l)/2) 2 . In figure 2, the cuvres of Eo n agaisnt a reveals that the 
Dirac-Hulthen spectrums are more evenly distributed with respect to the spatial dimensions; 
moreover, energy levels at higher exited states, tend to be bounded within [-1, 0]. 

Moreover, for large a, the energy eigenvalues (30) and (33) becomes imaginary. In 
particular if a > 2n+ 4 D _ 3 , the Klein-Gordon energy (30) becomes imaginary for the state 
E Kn (n, a, D). Similarly, if a > 2 (l + + p 2 ) /p 2 , where p 2 = (n + (D - l)/2) 2 , then the 
Dirac energy (33) becomes imaginary for the state Ec, n (n, a, D). 

4.2 Case 2 : Dependence of E n on D 

To investigate the dependence of the relativistic Hulthen energies on the dimension, we follow 
a recent study [40, 41] where the property of continuous dimension is used in studing the 
bound states of a quantum system with special potential. Thus, we assume the dimension D 
is continuous and plot some curves of E n against a as shown in figures 3 and 4 for the cases 
of Klein-Gordon and Dirac spectrum respectively. 

It is interesting to note that the set of points in figure 3 indicates that the Klein-Gordon- 
Hulthen energy preserves its discrete nature regardless of the the nature of D and the value of 
a. This is evident from the fact that in Eq. (32), if a is sufficiently small and D is sufficiently 
large, the energy Ek„ tends to -4, thus becoming independent of D for any given value 
of n. On the other hand, figure 4 reveals that the Dirac-Hulthen spectrum decreases and 
converges to a point with increase in the dimension and also for large a, there is a posssiblity 
of degenerate energy with respect to n. 



5 Concluding remarks 

We have studied the D-dimensional Dirac equation with the Hulthen potential within the 
frame work of an exponential approximation of the centrifugal term and position-dependent 
mass. The approximate eigenvalues obtained were found to be in good agreement with pre- 
vious works [17]; the eigenfunctions are also obtained in form of the Jacobi polynomials. The 
dependence of the energy spectrum on the screening parameter a and dimensions D were 
studied for both the Z)-dimensional Klein-Gordon- and the Dirac-equations with the Hulthen 
potential. We note that similar to the non-relativistic case, the Klein-Gordon-Hulthen energy 
levels for some adjacent dimensions intersect at some values of a, thus resulting into a de- 
generacy of energy eigenvalues with respect to the dimensions. However, the Dirac-Hulthen 
spectrums are more evenly distributed with respect to a and D. Finally, we observed that 
because of the non- linear terms in the energy equations (32) and (33), not all values of a yield 
a real eigenvalues for a given quantum state. 
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FIGURES 
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Figure 2: The curves of the Dirac- Hulthen energy spectrum against a for some dimensions 
and some exited states. 
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Figure 3: The curves of the Klein-Gordon-Hulthen energy spectrum against D for some exited 
states. The points '+', V and 'x' denote n = 3, 4 and 5 repectively 
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Figure 4: The curves of the Dirac-Hulthen energy spectrum against D for some exited states. 
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